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We present an in-depth analysis of the atomic and electronic structure of the quasi one-dimensional 
(ID) surface reconstruction of Ga on Si(112) based on Scanning Tunneling Microscopy and Spec- 
troscopy (STM and STS), Rutherford Backscattering Spectrometry (RBS) and Density Functional 
Theory (DFT) calculations. A new structural model of the Si(112)6 x 1-Ga surface is inferred. It 
consists of Ga zig-zag chains that are intersected by quasi-periodic vacancy lines or misfit disloca- 
tions. The experimentally observed meandering of the vacancy lines is caused by the co-existence of 
competing 6x1 and 5x1 unit cells and by the orientational disorder of symmetry breaking Si-Ga 
dimers inside the vacancy lines. The Ga atoms are fully coordinated, and the surface is chemically 
passivated. STS data reveal a semiconducting surface and show excellent agreement with calculated 
Local Density of States (LDOS) and STS curves. The energy gain obtained by fully passivating 
the surface calls the idea of step-edge decoration as a viable growth method toward ID metallic 
structures into question. 

PACS numbers: 68.35.-p, 68.37.Ef, 73.20.At, 81.07.Vb 



I. INTRODUCTION 



Nature only provides few one-dimensional (ID) elec- 
tronic systems, such as carbon nanotubes, 1 organic 
charge transfer salts, and inorganic blue bronzes (see 
for example the discussion in Ref. 0). Electrons con- 
fined to one dimension are fundamentally different from 
the quasi-particles of Fermi liquid theory^ In ID, even 
in the case of arbitrary low interaction strength, the 
single-particle description of the system breaks down and 
must be replaced by a description based on collective 
excitations?^ Experimental realization and verification of 
this Luttinger liquid phenomenon continues to capture 
the imagination of physicists, especially since the fab- 
rication of structurally uniform ID nanostructures now 
appears to be within the realm of possibilities. 

A very intuitive approach to produce ID systems is 
to utilize high index silicon surfaces. 4 Based on the con- 
cept of metal-adatom step-edge decoration, deposition of 
a submonolayer amount of metal atoms onto a stepped Si 
surface is expected to result in a single domain of quasi 
ID, metallic atomic wires, i.e. an atom wire array^ In 
contrast to for example carbon nanotubes, such a sin- 
gle domain surface quantum wire array would be eas- 
ily accessible to both nanoscopic and macroscopic tech- 
niques such as Scanning Tunneling Microscopy and Spec- 
troscopy (STM and STS), photoemission spectroscopy, 
and (surface) transport measurements. In addition, the 
coupling strength between the atom wires can be tuned 
by changing the miscut angle of the vicinal Si surface, 



i.e. adjusting the separation between the wires£ 

Indeed such single domain ID metallic systems have 
been produced on high index Si surfaces ^Ll&l How- 
ever, the metal adatoms in these studies are generally 
not adsorbed at the step edges. Instead, rather com- 
plicated reconstructions are formed with chains of metal 
atoms that are incorporated into the (lll)-like terraces 
in the unit cell^iSiii This questions the idea of form- 
ing atom wires via step-edge decoration on Si. Further- 
more, it was noticed^ that these ID atomic-scale sys- 
tems all exhibit intrinsic spatial disorder in the atomic 
structure, which will have important consequences for 
electronic transport in these systems. It should be noted 
that in the case that step-edge decoration does not occur 
in these studies, the miscut or vicinal orientation of the 
Si surfaces mainly serves to create a single domain sur- 
face reconstruction; similar or "parent" reconstructions 
exist on the corresponding planar surfaces. These single 
domain quantum wire arrays have been studied success- 
fully with Angle Resolved Photoelectron Spectroscopy 
(ARPES) 5 ! 9 ! 13 : 14 ! 15 and transport measurements^ 

In this paper, we investigate the formation of Ga chains 
on the vicinal Si(112) surface. A structural model for this 
interface was devised by Jung, Kaplan and Prokes (the 
JKP-model)i 17 i 18 i 19 i 20 i 21 i 22 The unit cell of the bulk ter- 
minated vicinal Si(112) surface contains a double-width 
(lll)-like terrace with single (lll)-like steps. Based 
on Low Energy Electron Diffraction (LEED) and Auger 
Electron Spectroscopy (AES) experimentsjiii&i 9 " it was 
proposed that Ga atoms adsorb at the step edges of the 
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bulk terminated unit cell, thus forming atom rows along 
the [110] direction. Missing Ga atoms or vacancies in 
these Ga rows align into quasi ID vacancy lines that 
run orthogonal to the Ga rows, resulting in the observed 
6x1 periodicity. In this model the Ga coverage is g of a 
Si(lll) bilayer, or five atoms per 6x1 unit cellii^ Later 
STM experiments by Baski et a/i22*2i seemed to confirm 
this model. These authors observed a well-ordered ar- 
ray of single-atom rows with a regular row spacing 9.4 A, 
equal to the step-edge spacing of bulk terminated Si(112); 
see for example Fig. 1 of Ref. |2(1 A side view and a top 
view representation of this JKP-model is shown in Fig. ^ 




FIG. 1: (a) Topview and (6) side view of the JKP-model 
of the Si(112)6 x 1-Ga surface. In (a) the 6x1 unit cell is 
indicated dotted. In (6) a (111) plane is indicated. Si atoms: 
light, Ga atoms: dark. 

As a consequence of the three-fold coordination of the 
adsorption sites, the trivalent Ga atoms are fully coordi- 
nated. There are no partially filled dangling bonds on the 
Ga atoms and the covalently bonded Ga atoms would not 
contribute any state density near the Fermi level. How- 
ever, an interesting feature which has remained largely 
unnoticed in literature is the fact that within the JKP 
model, there should exist a metallic dangling bond wire 
that is located on the row of Si surface atoms located 
in between the Ga rows. However, the predicted ID 
metallicity turned out to be unstable with respect to a 
Jahn- Teller distortion, leaving only one unpaired electron 
per 6x1 unit cell. Interestingly, the resulting electronic 
structure implied the existence of conduction channels 
orthogonal to the Ga chains^ 

We have performed a detailed Scanning Tunneling Mi- 
croscopy (STM) study of the Si(112)6 x 1-Ga surface. 
Because of the unprecedented resolution in the STM im- 
ages of the Si(l 12)6 x 1-Ga surface, a detailed investi- 
gation of the atomic structure of the Si(112)6 x 1-Ga 
surface could be carried out. Extensive Density Func- 
tional Theory (DFT) calculations have been performed 
to explore new candidate structural models. Theoretical 
STM images were calculated for the new structures and 
compared with the experimental STM images. From a 



detailed analysis of all experimental and theoretical in- 
formation, a new structural model for the Si(112)6x 1-Ga 
surface emerged, which shows excellent agreement with 
the experimental evidence. It contains two Ga atom rows 
amounting to a total of ten Ga atoms per 6x1 unit 
cell, consistent with RBS experiments. The two Ga rows 
form zig-zag chains while quasi-periodic vacancy lines in- 
tersect these Ga chains. The observed meandering of 
the vacancy lines can also be fully explained within this 
model. STS measurements show that the surface is semi- 
conducting, and are consistent with our band structure 
calculations and theoretical STS simulations. This paper 
presents a follow-up of an initial report^ with new data 
and provides a more detailed and in depth analysis, in- 
cluding a detailed comparison between spatially resolved 
STS and theoretical local density of states (LDOS) cal- 
culations. 



II. EXPERIMENTAL AND THEORETICAL 
PROCEDURES 

Experiments were carried out in an ultra-high vacuum 
system with a base pressure < 5 x 1CP 11 mbar. The sys- 
tem was equipped with a Ga effusion cell, direct current 
sample heating facilities, an Omicron variable tempera- 
ture STM and a LEED system. An n-type Si(112) wafer 
(~ 5 x 10 14 cm -3 , orientation ± 2° of the nominal (112) 
orientation) was cut into (10x2)-mm 2 samples and rinsed 
in acetone and isopropanol. After introduction into UHV 
the samples were degassed at 775 K overnight and sub- 
sequently the sample temperature was slowly raised to 
1025 K and kept there for 4 hrs. Next, the sample was 
flashed at 1475 K to remove the native oxide. During 
resistive heating, the current was directed parallel to the 
nano-facets of the clean (112) surface (i.e. in the [110] 
direction) in order to avoid current-induced step bunch- 
ing. The surface reconstruction was prepared in two dif- 
ferent ways. In the "one-step" procedure, Ga was de- 
posited with the Si substrate held at 825 ± 50 Ki24 In 
the "two-step" procedure, Ga is deposited onto a Si(112) 
substrate kept at room temperature. After Ga deposition 
the sample was annealed at about 825 ± 50 K to form 
the 6x1 reconstruction and to desorb excess Ga atomsiSS 
Both surface preparation procedures resulted in identical 
LEED patterns and STM images. The pressure remained 
below 2xl0 -10 mbar during sample preparation. The 
sample temperature during sample preparation was mea- 
sured using an optical pyrometer. STM and STS exper- 
iments were performed at room temperature and at low 
temperature (~ 40 K) using etched tungsten tips. STM 
images of the filled and empty electronic states were ob- 
tained with a constant current between 0.05 and 0.2 nA 
and bias voltages between 1 and 2 V. STS data were 
acquired with a setpoint of 0.3 nA at 1 V. 

RBS experiments were carried out at the AMOLF in- 
stitute in Amsterdam to determine the amount of Ga 
atoms per surface unit cell. A normal incident 2.0 MeV 
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He + ion beam from a Van de Graaf accelerator was 
backscattered from the Si crystal and detected at a 
backscattering angle of 165°. The beam current was typ- 
ically about 20 nA. 

The atomic and electronic structure of new candidate 
structural models for the Si(112)6 x 1-Ga surface, corre- 
sponding to Ga-coverages ranging from 5 to 11 Ga atoms 
per 6x1 unit cell were explored using an efficient local- 
orbital (LO) DFT technique (the Fireball96 code)^ 
In these calculations, we have used a minimal atomic- 
like basis set using the following cut-off radii (R c ) for the 
definition of the Fireball96 orbitalsS&: R c (Si) = 5.0, 
R c (Ga) — 5.2. For the most promising structures, Plane- 
Waves (PW) DFT calculations (Castep code)Sl were 
also performed to check the validity of the Fireball96 
findings. In these PW calculations, we have used a ki- 
netic energy cut-off E c of 200 eV for the definition of 
the PW basis set, and 4 special fc-points for the Brillouin 
zone sampling (test calculations with 250 eV and 8 spe- 
cial fc-points were also performed). In both the LO and 
PW calculations we have used a slab of 11 Si layers with 
hydrogen atoms saturating the bonds of the deeper Si 
layers (see Fig. QJ. 

Using the DFT local-orbital hamiltonian of the sur- 
face together with non-equilibrium Keldysh Green func- 
tion techniques we calculated theoretical STM im- 
ages for the new relaxed atomic structures. The theoret- 
ical images were then compared with the experimental 
STM images. In our approach, we divide the total hamil- 
tonian^ H, of our tip-sample system into three parts, 
H = H t + H s + H int , H t , H s and H int referring to the 
tip, sample and their interaction. H s is obtained from the 
FiREBALL-code used to calculate the Si(112)6 x 1-Ga sur- 
face; H t is calculated using the same DFT local-orbital 
code for a W-tip having a pyramid with four atoms, at- 
tached to a W-(100) surface; H int is obtained using a 
dimer approximation, whereby the different tip-sample 
hopping interactions, T ts , are calculated from the dimer 
formed by the respective tip and sample atoms whose in- 
teraction we want to obtain (it is shown in Ref. |30j that 
this approximation yields a good description of the STM 
images if orbitals with long-range tails are used in the 
hopping calculations). A more detailed description of our 
procedure to obtain theoretical STM images can be found 
in Refs. I30ll3ll Making use of the total Hamiltonian and 
the Keldysh Green-function techniques, we can calculate 
the tunneling current from the following equation 28 

Aire f°° » 
1= — j dLoTr[T ts p^)DUu)T st p u (u)D a tt {u)\ 

x(/ t (w)-/ a (w)) (1) 

where 

^[i-ik^HT^LMr 1 (2) 

and 

bl = [i-f ta rMf st g a tMr 1 (3) 



include all the interface multiple scattering processes. 
Tr stands for the Trace of the current matrix, g^ r ^ and 
are the advanced (retarded) Green-functions of the 
sample and the tip, respectively (calculated taking T ts — 
0); p ss and ptt are the sample and tip density of states 
(also for T ts — 0); and ft (f 8 ) the Fermi distribution 
functions. 

In the tunneling regime, T ts is very small and D r ss 
and £>tt can be replaced by /. In this limit, for zero- 
temperature, we recover the following equation: 

I = —r-l dLuTi[T ts p ss (cu)T st ptt(uj)] (4) 
11 Je f 

which we have used to calculate the STM images of 
the different surface structures. 

We should comment that the detailed comparison be- 
tween theory and experimental results that we intend 
in this work requires the use of equation J3J instead of 
other simpler approaches (like the Tersoff-Hamann for- 
malism) that are common in the literature. Our method 
includes a realistic description of the geometry and the 
full electronic structure of the tip, and incorporates quan- 
titatively the influence of the tunneling parameters (bias 
and current conditions) and the tip-sample distance^S^i 
This quantitative accuracy, crucial to understand the 
contradictory experimental results in terms of contrast 
and symmetry of the STM images of an apparently sim- 
ple system like 0/Pd(lll)-2x2;2iis necessary in our case 
to discriminate among all the different surface structures 
that have been analyzed in this work. Notice, in particu- 
lar, that we show below that our proposed model is fully 
compatible with the STM images by Baski el alm& pro- 
vided that their tunneling parameters are used in the sim- 
ulation of the STM images. On top of these advantages, 
we have to mention that our approach does not require 
a significantly larger computational time than other sim- 
pler methods, as equation J3J provides a very compact 
procedure for calculating the tip-sample tunneling cur- 
rent that takes full advantage of the LDOS (p" ss and p" tt ) 
obtained from our DFT calculations. 



III. STM OBSERVATIONS 

The high index Si(l 12) surface is tilted 19.5° away from 
the (111) surface towards (001). But the pristine Si(112) 
surface is not thermodynamically stable, and breaks up 
into approximately 10 nm wide nano-facets of recon- 
structed (111)- and (337)-like planesi^S An STM image 
of pristine Si(l 12) is shown in Fig.^a). One might ex- 
pect that metal deposition on this surface would result in 
the formation of metallic nanowires in these prepatterned 
grooves. However, it was shown by Baski el alm& that 
upon deposition and post-annealing of a sub-monolayer 
amount of Ga the faceted Si(112) surface undergoes a 
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massive restructuring. It returns to its basal (112) orien- 
tation, reconstructing as described in the introduction. 
This preparation procedure of deposition and postan- 
nealing of the surface resulted in a reproducible self- 
limiting surface reconstruction with a 6 x 1 unit cell^ 
A large scale STM image of Ga covered Si(112) is shown 
in Fig. The nanoscale facets have developed into 

large anisotropic (112) terraces that can extend for up 
to microns along [110] direction but are less than 100 
nm wide. On the terraces, the vacancy lines appear as 
dark trenches which run perpendicular to the step edges 
present in this image. Closer inspection reveals that these 
vacancy lines are not exactly straight, but their position 
fluctuates around an average position. As discussed be- 
low, this is due to the coexistence of 6 x 1 and 5x1 
units in the surface and due to the presence of intrinsic 
fluctuations in the vacancy finest 

Detailed atomic resolution STM images were acquired 
to investigate the atomic structure of the Si(112)6 x 1- 
Ga surface. In Fig. |2| an atomic resolution empty state 
STM image is shown. Note that this particular area of 
the surface shows both 5x1 and 6x1 unit cells as indi- 
cated in the figure. Two parallel atom rows are observed 
per unit cell, running in the [110] direction, intersected 
by the quasi-periodical vacancy lines. Comparing this 
image with the results of Baski et aZ.jSSl we observe the 
same spacing of the brightest atom rows (i.e. 9.4 A). Fur- 
thermore, the mixed periodicities and the similar LEED 
pattern (see Refs. I20I and Il8l respectively), indicate that 
the same surface reconstruction is studied here. Con- 
sequently we conclude that the brightest atom rows in 
Fig. are the same atom rows as imaged by Baski et 
al~ (henceforth, the "step-edge Ga row"). But in ad- 
dition we observe a 2 nd atom row lying in between the 
brighter rows. In terms of the JKP-model, this row of 
atoms could be interpreted as the Si dangling bond row 
which might form a quasi ID band. However, these two 
parallel atomic lines clearly form a zig-zag pattern as in- 
dicated in figures |21 (see also Fig. IT3|) , which results in 
a structural asymmetry in the vacancy line. This is in 
contradiction with the JKP-model, which implies mirror 
plane symmetry with respect to the (110) plane in the 
vacancy line. 

In Fig. 0] a set of registry aligned dual bias images is 
presented. These images have been recorded simultane- 
ously on the same area of the surface, but with opposite 
tunneling bias polarities resulting in a set of spatially cor- 
related empty and filled state images. In this case, the 
empty state image, Fig. Ufa), has suffered from a slight 
decrease in resolution, as compared with Fig. |3 but the 
asymmetry in the vacancy line is still visible. In the 
filled state image, Fig. Hffe), a relatively big, symmetric 
protrusion prevents a detailed observation of the atomic 
structure in the vacancy line. As in the empty state im- 
age, two parallel atom rows are also visible in the filled 
state image. They form a ladder structure instead of the 
zig-zag pattern of the rows observed in the empty state 
image. 




FIG. 2: (a) STM image of pristine Si(112). (6) STM image 
of the Ga covered Si(112) surface. Tunneling conditions: 1.5 
V, 0.1 nA and 2 V, 0.1 nA, for (a) and (6), respectively. 



RBS measurements were performed to determine ex- 
perimentally the amount of Ga at the surface. Integra- 
tion of the Ga peak in the backscattered He spectrum 
yielded an amount of 9 ± 1 Ga atoms per 6x1 unit cell, 
to be compared with 5 Ga per 6x1 unit cell for the 
JKP-model of Fig. [TJ 

In summary, these experimental results (STM and 
RBS) consistently show that the step-edge decorated 
JKP-model of the Si(112)6 x 1-Ga surface is at variance 
with the new experimental observations. Consequently, 
the intuitive idea of metal adatom step-edge decoration 
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FIG. 3: Empty state STM image of the Si(112)n x 1-Ga sur- 
face. In this particular surface area both 5x1 and 6x1 unit 
cells are present, as indicated. Atomic positions are indicated 
with black dots. The position of the vacancies in both atomic 
rows is indicated with white dots. Tunneling conditions: 1.5 
V, 0.2 nA. 




FIG. 4: (a) Empty state and (6) filled state dual bias STM 
image of the Si(112)6 x 1-Ga surface. In (6) the ladder struc- 
ture is indicated. Tunneling conditions: ± 1 V, 0.05 nA. The 
inset shows a filled state image with a slightly lower resolu- 
tion. Tunneling conditions: -2 V, 0.1 nA. 

does not seem applicable for the Ga/Si(112) interface. 

IV. STM IMAGE SIMULATIONS 

Extensive DFT calculations were performed to identify 
the precise atomic structure of the Si(112)6 x 1-Ga sur- 



face. Using the DFT local-orbital hamiltonian of the sur- 
face together with non-equilibrium Keldysh Green func- 
tion techniques ^LStL we calculated theoretical STM im- 
ages of these most promising structures, which were then 
compared with the high resolution experimental STM im- 
ages. The different Si(112)6 x 1-Ga structures analyzed in 
this paper, with Ga-coverages ranging from 5 to 11 Ga 
atoms per 6x1 unit cell, have been generated starting 
from the JKP-model (see Fig. in the following way: 
(a) replacing some of the Si atoms in the Si-dangling- 
bond row by Ga atoms (hereafter referred to as Ga ter- 
race atoms); (b) replacing some of the Ga atoms at the 
step edge by Si atoms; (c) considering also the replace- 
ment of Si or Ga atoms on the step-edge and terrace 
rows by vacancies and the addition of Ga or Si atoms in 
the vacancy lines. In total, more than 40 new structures 
were fully relaxed, their surface energies and electronic 
structures calculated, and their corresponding theoreti- 
cal STM images obtained. In the following, we compare 
the theoretical STM images of structures with the lowest 
total energies; a detailed chemical potential analysis of 
the total energies is deferred to Section IVT1 

Fig. [SJa) and (b) show two examples of simulated 
STM images for some of these structural models, with a 
top view of the corresponding atomic structure superim- 
posed. Fig. a) corresponds to a structural model that 
contains, per 6x1 unit-cell, 6 Ga atoms in the terrace 
row, and 5 Ga atoms plus a vacancy in the step-edge row; 
in Fig.|3J&) there are 5 Ga and one Si in the terrace row, 
and 5 Ga plus a vacancy in the step-edge row. The simu- 
lated STM images for the different structural models are 
compared in detail with the experimental high-resolution 
STM images. For example, Fig.[SJa) (filled state) is sim- 
ilar to the filled-state STM image obtained in Ref. l2(i 
also the empty-state image of Fig- Et is m good agree- 
ment with the experimental STM image shown in Fig. [3] 
However, a detailed analysis of both empty and filled 
states images, as well as registry aligned dual images (e.g. 
Fig. |3J , reveal that these models present some inconsis- 
tency with the experimental high-resolution information. 
For example, the structural model of Fig. a) is sym- 
metric with respect to the vacancy line, in disagreement 
with Fig. |3 in the case of Fig.[Sf &) registry aligned dual 
bias STM images show that the bright protrusion in the 
filled state image is located in the vacancy line, aligned 
with the brighter Ga row of the empty state image, while 
in the simulated filled state-image it appears in between 
the two Ga rows. 

Thus, a detailed comparison of the theoretical STM 
images for the different structural models with the ex- 
perimental STM images was performed. From this anal- 
ysis, we concluded that the correct atomic model for the 
Si(112)6 x 1-Ga is the one shown in Fig. [HJ In this new 
structural model there are 10 Ga atoms per 6x1 unit- 
cell (to be compared with the RBS determination of 9 ± 1 
Ga atoms), forming two parallel rows, in a zig-zag con- 
figuration (see also Fig.|3|). The upper row of step-edge 
Ga atoms adsorbed at the (lll)-like step is equivalent to 
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FIG. 5: Simulated empty (left) and filled (right) state STM 
images of some of the structural models analyzed, with a ball 
and stick representation (top view) superimposed on top of 
the STM images, (a) structural model with, per 6x1 unit- 
cell, 6 Ga atoms in the terrace row and 5 Ga atoms in the 
step-edge row; (b) a structural model with, per 6x1 unit-cell, 
5 Ga and 1 Si in the terrace row, and 5 Ga in the step-edge 
row. Ga atoms: dark, Si atoms: light. 

the Ga row in the JKP-model. But the Si dangling bond 
row in the JKP-model has been replaced by a second row 
of Ga atoms (henceforth, the "terrace Ga row"). Each 
Ga-row contains 5 Ga atoms per 6x1 unit-cell, i.e. there 
is a Ga-vacancy in each row. These vacancies are placed 
at adjacent sites in the (zig-zag) two rows, giving rise to 
an asymmetry in the vacancy line, see Figs. and [21 

Inside the vacancy lines, missing Ga atoms expose the 
underlying Si atoms. These Si atoms rebond forming Si- 
Si dimers on the terraces and Si-Ga dimers along the step 
edges in each unit cell. Specifically, by rotating a step- 
edge Si atom toward the step-edge Ga row, this Si atom 
can rebond to two neighboring Si atoms and a step-edge 
Ga atom, forming a Si-Ga dimer with the latter. The 



rebonding of the Si atoms in the vacancy line also im- 
plies that the Ga vacancies on both Ga rows must be 
aligned. In contrast with the tetravalent Si atoms in the 
step-edge decorated JKP-model, both the trivalent Ga 
atoms in the three-fold adsorption sites on the (lll)-like 
terraces, and the exposed Si atoms inside the vacancy 
lines present no unsaturated dangling bonds; the result- 
ing structure (Fig. EJ is fully passivated and the surface 
is semiconducting. 




FIG. 6: Ball and stick representation of the energy minimized 
structure for the Si(112)6 x 1-Ga surface; topview (a), and 
sideview (b). In (a) a 6 x 1 unit cell is indicated and in (b) a 
(111) plane is indicated. Si atoms: light, Ga atoms: dark. 

The calculated theoretical STM images corresponding 
to this zig-zag model are shown in Fig.[7| with a top view 
of the structural model superimposed on top. Both the 
empty state and filled state images are in excellent agree- 
ment with the experimental ones (see Fig. [3J|. It shows 
that the two atom rows imaged in the empty state STM 
images, are indeed the step-edge Ga row and the terrace 
Ga row, ruling out the formation of a Ga-atom step-edge 
decorated structure. In addition, the asymmetry in the 
vacancy lines observed experimentally in the empty state, 
is neatly reproduced in the simulated STM images. In the 
simulated filled state image, fuzzy lines with a big, sym- 
metric protrusion inside the vacancy line are observed, 
in agreement with the experimental images. Fig. 
clearly shows that the big protrusion corresponds to the 
Ga-Si dimer. Furthermore, it shows that the fuzzy lines 
are originating from a Si-Ga bond on the (lll)-like ter- 
race. They form a ladder configuration, in agreement 




FIG. 7: Simulated empty (a) and filled (b) state STM images 
of the zig-zag model of Fig. El A top view of the ball and stick 
representation is superimposed on top of the STM images. 
Ga atoms: dark, Si atoms: light. Zig-zag symmetry and 
ladder symmetry indicated with white bars in (a) and (b), 
respectively. Tunneling bias 2 V (a), and - 1.3 V (b). (c) 
Side view of the proposed model. The (111) plane is indicated 
with a dotted line, (d) Simulated filled state image, -2 V. 

with the atomic resolution experimental image in Fig. 0] 
The only feature which was not reproduced is the slightly 
higher apparent height, in the empty state experimental 
images, of the Ga atoms in the two terrace Ga rows di- 
rectly adjacent to the vacancy line. Finally, we mention 
that changing the tunneling conditions in the simulated 
STM images (tip-sample distance, voltage), the experi- 
mental STM images of Ref.|2(]can be recovered, as shown 
in Fig.Etd). 

V. SPECTROSCOPY 

We also have studied this surface reconstruction with 
scanning tunneling spectroscopy. While imaging the sur- 
face with constant tunneling current, at every third data 
point an / — V curve is measured with the feedback loop 
switched off during this I — V measurement. In Fig. we 
have averaged I — V curves measured on the upper atom 
rows, the lower atom rows, and on the vacancy lines sep- 
arately (the respective areas being determined from the 
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FIG. 8: Semilogarithmic (a) and normalized derivative (6) 
plots of I — V-curves averaged over the step-edge (dotted) 
and terrace (dashed) atom rows, and the vacancy lines (solid), 
respectively. STS setpoint: 1 V, 0.3 nA. The inset in (a) 
shows the bandstructure inferred from the data. 

empty state STM image). This results in three curves, 
representing the electronic structure on the terrace and 
step-edge Ga rows, and the electronic structure inside 
the vacancy line. At bias voltages below the bulk con- 
duction band minimum, the tunneling current is limited 
by thermionic emission, as is evident from the linear in- 
crease of the log(I) — V curve, 33 up to the conduction 
band minimum (CBM) at 0.7 V, see Fig. ^a). Conse- 
quently, the bulk valence band maximum (VBM) at the 
surface should be located at ~ —0.4 V, implying an up- 
ward band bending of ~ 0.5 eV for this n-type specimen 
(10 15 cm~ 3 ), E f - E V bm and E CBM - E f being - 0.4 
and ~ 0.7 eV at the surface, respectively, as shown in 
the inset of Fig. Ufa). These data are consistent with 
the measured surface photovoltage in Refs. 12^1341 In 
Fig. ffife) we have plotted the normalized derivative of 
the three I — V curves (i.e. g^ v ), originating from the 
three different areas within the unit cell. These tunnel- 
ing spectra are proportional to the local density of states 
(LDOS) at the respective areas of the samp le surface over 
which the averaging took place (see Refs. I35l36l37l) . In 
all curves, there is no DOS at the Fermi level, but a 
gap of exists between the filled and empty state bands 
showing that indeed the surface is semiconducting. The 
tunneling spectra on the two Ga atom rows appears to be 
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FIG. 9: (a) |£-map of measured STS curves at -0.31 V. STS 
setpoint: 1 V, 0.3 nA. (6) Corresponding topographic STM 
image. Tunneling conditions 1 V, 0.3 nA. 



similar in shape, whereas the tunneling spectra inside the 
vacancy line deviates from the former two. The leading 
edge of the total tunneling spectrum in the filled state 
spectrum arises from a state (B) that is located mainly 
inside the vacancy lines at about ~ —0.3 eV. The two 
Ga rows feature a broad filled state at higher binding en- 
ergy (C). In the empty state tunneling spectra the two 
atom rows show a small shoulder (A) at about ~ 0.6 
eV, just below the bulk CBM. Thus a surface band gap 
of ~ 0.9 eV is deduced. These experimental data are 
entirely consistent with the presence of two equivalent, 
threefold coordinated rows of Ga atoms at these posi- 
tions, fully passivating the surface. Notice that we in- 
ferred the surface band gap from the separation between 
peaks A and B using the peak position or centroids, and 
not the onsets. The justification for this procedure comes 
from a detailed comparison with theoretical STS data, as 
will be discussed below. 

In addition, we have constructed a ^-map of the STS 
measurements. In a J^-map, the value of the derivative 
of the I — V-curves at a certain voltage V is plotted as a 
two-dimensional image, with the x and y coordinates cor- 
responding to the topographic STM image. In Fig.|§fa) 
the derivative of the I — ^-curves at -0.31 V is plotted, 
the corresponding empty state STM image is shown in 
Fig. Indeed, the largest slope in the I— ^-curves at 
-0.31 V is located inside the vacancy lines (i.e. here the 



largest increase in tunneling current is observed, corre- 
sponding to the largest LDOS as compared to the LDOS 
at this specific energy at other locations on the surface). 
Only a very small intensity variation is observed perpen- 
dicular to the atom rows, consistent with the similar filled 
state tunneling spectra on the two Ga rows in Fig.ffia,). 

From the structural model and its spatially re- 
solved DOS, theoretical STS curves were calculated, see 
Fig. HOf b). As for the STM images, these results are 
obtained using the LO-DFT Fireball96 hamiltonian of 
the surface and the Keldysh Green function approach. 
The corresponding calculated LDOS of the zig-zag model, 
averaged over different areas is shown in Fig. HOf a") (a 
broadening of 0.1 eV has been used). The calculated STS 
curves shown in Fig. HOf 6) were obtained by placing the 
tip over the respective areas, calculating the current as 
a function of a voltage sweep using the calculated LDOS 
(Fig. 1101) . and averaging over the areas of interest. The 
Fireball96 local orbital calculations employ a minimal 
basis set, resulting typically in band gaps that are too 
large. Nonetheless, excellent qualitative agreement ex- 
ists between the calculated STS curves and normalized 
derivatives of the experimental STS curves (Fig. Et^))- 
The calculated LDOS and STS curves confirm that the 
large peak B just beneath the band gap indeed is mainly 
associated with states that are located on the Si-Ga dimcr 
inside the vacancy lines. The two Ga rows contribute 
almost equally to a broad peak in the DOS at higher 
binding energy (C) and a small shoulder in the DOS just 
above the band gap (A), in full agreement with the nor- 
malized derivative of the tunneling spectra in Fig. 
The empty state DOS is very similar for both Ga rows. 
Consequently the ~ 0.4 Ahigher appearance of the Ga 
atoms at the step edge in the empty state image is due 
to their on average higher atomic positions, and thus the 
empty state STM image reflects the real surface topog- 
raphy at these voltages. 

We have calculated the surface bandstructure using 
both the LO and PW-DFT methods, within the LDA for 
exchange-correlation contributions. While the LO calcu- 
lation overestimates the value of the bulk band gap, the 
PW calculation typically underestimates the band gap. 
In the LO bandstructure (not shown) a surface state band 
gap of 1.2 eV is obtained between the states A and B, 
while the separation between peaks A and B in the calcu- 
lated STS is close to 1.4 eV. This suggests that in order to 
measure the band gap, it is reasonably accurate to use the 
peak positions instead of the (poorly defined) peak onsets 
in the experimental Jy^y curves. The observed splitting 
between A and B in the experimental curves is ~ 0.9 eV 
indicating that the experimental band gap is ~ 0.8 eV. 
Fig. El shows the bandstructure as calculated with the 
PW-DFT code, showing a surface band gap of 0.77 eV 
between states A and B4£ The calculations place filled 
state B slightly above the VBM, in agreement with the 
experimental observation. On the other hand, the empty 
surface state A is located at or slightly above the CBM 
according to the PW-DFT calculations, while experimen- 
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tally state A appears slightly below the bulk CBM; see 
Fig|H». The PW gap of 0.77 eV is comparable to the 
experimental band gap of 0.8 eV. However, the precise 
location of state A in the calculations directly affects the 
value of the band gap. 
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FIG. 10: (a) Calculated LDOS (Fireball96), averaged over 
the step-edge (dotted), terrace (dashed) Ga-rows, and va- 
cancy line (solid), respectively. The LDOS at the vacancy 
line is calculated as the average density of states of the Ga-Si 
dimer and Si-Si dimer. (b) Normalized derivative of the I—V- 
curves shown in (a). A broadening of 0.1 eV has been applied. 
Note the larger energy scale as compared to Fig.|H]due to the 
overestimation of the gap in the local orbital calculation. 




FIG. 11: Calculated band structure (Castep) of the zig-zag 
structural model. Surface states labelled A and B are men- 
tioned in the text. The shaded area shows the bulk projected 
bandstructure of the Si substrate. The inset shows the surface 
Brillouin zone probed. 



VI. 



COMPETING STRUCTURES 



The discussion above shows that, to elucidate the pre- 
cise atomic structure of a complex surface like Si(112)6 x 
1-Ga from the comparison of theoretical and experimen- 
tal STM images, it is necessary to use high-resolution ex- 
perimental STM images, including registry aligned dual 
bias information and STS data, combined with state of 
the art theoretical STM simulations. These simulations 
were performed on the subset of possible structures that 
were deemed most realistic on the basis of total energy 
considerations. In this section, we explore the relative 
stability of the various structures, which gives a more 
physical basis to the proposed structural model. 

In general, the precise stoichiometry of the surface is 
not known, and thus the analysis of the relative stability 
of different structural models requires the calculation of 
the surface energy F as a function of the different chem- 
ical potentials. 



Despite the fact that the Ga/Si(112) overlayer ap- 
pears to be quasi two-dimensional in atomic structure, 
the electronic structure of this overlayer is quasi one- 
dimensional. The Ga-atom induced surface band A dis- 
perses around the AT-point minimum in the upper part 
of the band gap. The dispersion of this band near the Ap- 
point yields an effective mass of m* ~ 1.48 m e along the 
X — r direction, and m* ~ 0.15 m e along X — K . This 
indicates a quasi-one-dimensional dispersion. This Ga- 
band is initially empty but could perhaps be populated 
in a controllable way, using a biased gate electrode, or a 
heavily n-type doped substrate, making the Si(l 12)6 x 1- 
Ga surface a promising system for the experimental study 
of electron transport in one dimensional atomic wires. 



A. Chemical potential analysis 

For the analysis of the relative stabilities of the vari- 
ous structures, we need to calculate the surface energy 
F = E t ot - HGaN Ga - fJ-SiNsi, where E tot is the total 
energy per unit-cell, [iQa , /J-Si are the Ga and Si chemical 
potentials, and Nc a , N$i are the number of Ga and Si 
atoms in the unit-cell. For we use the total energy 
(per atom) of bulk-Si (i.e. the surface is in equilibrium 
with the substrate). The value of ^,Ga is not determined 
by the substrate, but it can be estimated analyzing the 
experimental conditions (see below). 

Fig.lF2lshows the surface energy F — E to t — ^Ga^Ga — 
[J-SiNsi as a function of jJLQa- In this figure we use the 
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structural model of Fig. as reference, and plot F for 
some of the most promising models, as calculated with 
the PW code (Castep). In order to estimate the value 
of UGa we have to analyze the experimental conditions of 
the Ga deposition. In the "one-step" process the (6 x 1)- 
phase is formed under a Ga flux from the effusion cell 
with the sample held at a temperature of T = 825 K. At 
this temperature the incoming flux of Ga atoms is bal- 
anced by a flux of Ga atoms desorbing from the surface, 
thus establishing a quasi-equilibrium. This allows us to 
relate the chemical potential [iQ a in the effusion cell with 
the chemical potential fiQa at the sample. The chemical 
potential in the effusion cell may be approximated by the 
total energy of bulk-Ga, fiGa(pulk), i.e. the Ga vapor in 
the effusion cell is in equilibrium with the solid. 41 

Considering also the equilibrium between the sample 
and the Ga vapor in contact with the sample, we can 
estimate the chemical potential at the sample 

fJ-Ga = nca(bulk) - k B Tln{—) 

Ps 

where p c is the Ga vapor pressure in the effusion cell and 
p s the Ga vapor pressure at the sample. Since the effusion 
cell flux is proportional to its vapour pressure, p c , times 
the cell aperture area, and the sample flux is also pro- 
portional to its corresponding vapour pressure, p s , times 
the sample area, we conclude that (p c /p s ) ~ 10 2 (10 3 ), 
and fi Ga = nca-buik - 0.32(0.48) eV. In Fig. [T^ we see 
that for this range of pica the structural model of Fig. H3 
presents the lowest surface energy F. This result strongly 
supports our conclusion that the structural model for the 
Si(112)6 x 1-Ga surface is the one depicted in Fig. 




in surface free energy. The results shown in Fig. 1121 sug- 
gest, however, that the JKP-model might be stabilized 
for very low UGa values. We should stress that our the- 
oretical analysis has been directed to search for surface 
atomic structures that could explain the experimental 
results (STM and RBS) for the Si(112)6 x 1-Ga surface, 
and thus surface structures with lower Ga coverages, that 
should be favored for low nca values have not been ana- 
lyzed as thoroughly as those with coverages close to 9-10 
Ga atoms/6 x 1 unit cell. Nevertheless, we may perform 
a simple analysis, comparing the surface energy of the 
step-edge decorated JKP-model, with the surface energy 
of a simple hypothetical surface: half the surface is cov- 
ered with the structure of Fig. (i. e. both step-edge and 
terrace Ga rows) while the other half consists of clean 
Si(112). Both the step-edge decorated JKP system and 
this hypothetical half-half case present the same Ga cov- 
erage and thus the same behavior of F as a function of 
MGa [i.e.. the same slope in Fig. I12JI . This comparison 
reveals that the hypothetical case is lower in energy (by 
~ 0.7 eV/(10 Ga atoms), for all fiQa values, showing 
that the step-edge decorated case is unlikely to be sta- 
bilized at lower Ga coverage, and phase separation into 
bare Si(112) and the Si(112)6 x 1-Ga zig-zag surface will 
occur instead. Note that the facetting of the unstable 
Si(112) surface, which was not accounted for in this cal- 
culation, would increase this energy difference, making 
phase separation even more favorable compared to the 
step-edge decorated JKP-model. 

Another possible scenario for obtaining a step-edge 
decorated Ga row would be to use the experimentally ob- 
served fully passivated Si(112)6 x 1-Ga surface (Fig. EJ) 
as starting point, and try to kinetically stabilize a 
metastable step-edge decorated structure by selectively 
desorbing the Ga atoms from the terraces. Experimen- 
tally, this might happen in the 'two-step' preparation pro- 
cedure, see Section|n] We have studied this possibility by 
calculating desorption energies of terrace and step-edge 
Ga atoms from the Si(112)6 x 1-Ga surface (Fig. EJ. In 
particular, we have considered removing the Ga atoms 
close to the vacancy line as well as the replacement of 
those Ga atoms by Si atoms. In both cases the desorp- 
tion energies are lower by ~ 0.7 eV for step-edge Ga 
atoms than for Ga atoms on the terraces. This result 
suggests that a metastable step-edge decorated structure 
likely can not be achieved by thermally desorbing the Ga 
terrace atoms. 



FIG. 12: Surface energy as a function of the Ga chemical 
potential (Castep). The Ga chemical potential is plotted 
relative to the chemical potential in bulk Ga. 5Ga (filled 
squares) is the step-edge decorated JKP-model 

Comparing our model with the step-edge decorated 
JKP-model, an important difference is that the new 
structural model (Fig.|3 presents no partially- filled dan- 
gling bonds as discussed above. The stability of the new 
model is related to the full passivation of the substrate, 
removing all dangling bonds, and the associated decrease 



B. Intrinsic structural disorder 

In the experimental STM images, the vacancy lines are 
not exactly straight, but some meandering is observed, as 
shown in Fig.^5](see also Figs.|2J&) and Fig.[3J). The new 
structural model is able to fully explain the experimen- 
tally observed meandering. 1 — It was proposed by Erwin 
et alJL that this meandering of the vacancy lines could 
be explained by the co-existence of 6 x 1 and 5x1 unit 
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cells on the surface. This results in occasional steps in 
the vacancy lines, equivalent to the observed meander- 
ing of the dimer-vacancy lines on the Ge covered Si(OOl) 
surfaces^ We have analyzed the stability of the zig-zag 
structural model as a function of the longitudinal period- 
icity. Fig. El shows the surface energies F of this model 
for different periodicities: 5 x 1, 6 x 1 and 7x1. The 
5x1 surface corresponds to 4 Ga atoms in each Ga-row 
between vacancy lines, while the 7x1 surface presents 
6 Ga atoms in each row between vacancy lines. For our 
estimated range of [iQ a the 6x1 surface presents the 
lowest F, while the 5 x 1 is only 0.1-0.2 eV higher, per 
6x1 unit-cell. This small energy difference should lead 
to the experimental observation of 5 x 1 unit cells. Indeed 
these unit cells are frequently observed in the experimen- 
tal images, see Fig. ^0 It thus appears that the pre- 
dictions from the ID Frenkel-Kontorova model regard- 
ing the vacancy-line spacing in the step-edge decorated 
structure^! also apply to the quasi 2D zig-zag arrange- 
ment of Ga atoms presented here. This conclusion is 
not very surprising as it was concluded in Ref. |^] that 
the strain induced by the size difference of Ga and Si 
completely dominates the energetics of the periodicities. 
Apparently this conclusion still holds when an extra Ga 
row is added. 




FIG. 13: Detailed STM image, showing the two contributions 
to the meandering of the vacancy lines. Straight lines are 
drawn through domains with unit cells of the same size. Oc- 
casional jumps in the lines are due to different unit cell sizes, 
as explained in the text. Dots are placed on the Ga atom of 
the Ga-Si dimer at the vacancy of the step-edge rows, showing 
the fluctuations due to the intrinsic disorder of the random 
orientation of the Ga-Si dimer. 

However, careful investigation of the experimental im- 
ages reveals that the meandering of the vacancy lines as 




FIG. 14: Free energy as a function of the Ga chemical po- 
tential for 5x1, 6x1 and 7x1. The Ga chemical potential is 
plotted relative to the chemical potential in bulk Ga. Note 
the difference in scale on the abcissa as compared to Fig. 1121 

observed in Fig. is not only due to the competing lon- 
gitudinal periodicities. Instead, for large sections of the 
surface, the terrace Ga rows are perfectly periodic in the 
xl direction with (n-1) Ga atoms per terrace Ga row 
in the nxl unit cell. But in these ordered sections the 
number of Ga atoms in the step-edge Ga rows appears 
to fluctuate between n-2 and n. The proposed structural 
model perfectly explains these intrinsic fluctuations (i.e. 
fluctuations within a nxl domain); they are related to 
the orientation of the Si-Ga dimer in the step-edge rows. 
The twofold symmetry of the 112 substrate in the [110] 
direction is broken by the Si-Ga dimers, resulting in two 
degenerate orientations of these dimers. The energy as- 
sociated with interchanging the atoms of a Si-Ga dimer 
has been calculated™ to be less than 10 meV per 12 x 1 
unit cell. This small energy difference explains the ap- 
pearance of frequent meandering in the aligned vacancies 
in the step-edge Ga rows, thus accounting for the major- 
ity of the fluctuations in the vacancy lines observed in 
the experimental images. The absence of these fluctua- 
tions in the filled state images (compare Figs. EI a ) and 
(6) and also Figs.UJa) and (&)) is the result of the fact 
that the bright protrusion in the vacancy line appears 
in the center of the Ga-Si dimer, making its appearance 
insensitive to the orientation of the Ga-Si dimer. 



VII. SUMMARY AND CONCLUSIONS 

The 6x1 reconstruction of Ga on vicinal Si(112) was 
studied with STM, STS, RBS, and extensive DFT cal- 
culations. High resolution STM experiments revealed an 
asymmetry in the vacancy lines of the Si(112)6 x 1-Ga 
surface that is inconsistent with the JKP-model of step- 
edge decoration. STS measurements also rule out for- 
mation of quasi ID metal wires while RBS experiments 
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indicated a Ga coverage twice as large as previously in- 
ferred from the JKP model. Extensive DFT calculations 
were used to analyze the relative stability of more than 
forty structures, taking the chemical potential of the Ga 
adsorbate into account. Theoretical STM images were 
calculated for the most promising structures and com- 
pared in detail with the experimental STM images. 

A new structure emerged containing 10 Ga atoms per 
6x1 unit cell. The Ga atoms decorate the step edge and 
passivate the terrace atoms, thereby forming a zig-zag 
pattern. Excellent agreement between experimental and 
theoretical STM and STS data confirmed the validity of 
the proposed zig-zag model and demonstrate the power of 
such a comparison. Ga atoms are threefold coordinated 
and Si dangling bonds are all passivated so the surface 
is semiconducting. The "broken bond orbitals" inside 
the vacancy lines rebond to form Si-Ga and Si-Si dimers. 
The observed meandering of the vacancy lines originates 
from thermal fluctuations between the two symmetry- 
degenerate orientations of the Si-Ga dimer, in conjunc- 
tion with thermal fluctuations between competing 6x1 
and 5x1 units. 

While step-edge decoration of vicinal metal surfaces 
works*^ the observed drive toward chemical passivation 
suggests that step-edge decoration of vicinal semiconduc- 
tors is not a viable method to produce ID metal wires. 
Although this general conclusion remains to be tested 
further, it is clear that partially-filled dangling bonds on 
the terraces of vicinal surfaces are always greatly reduced 
in number or eliminated altogether in the reconstruction. 



As shown in this paper, predictive calculations along 
these lines should always take into account the chemi- 
cal potential of the adsorbate. The latter depends on 
the experimental preparation conditions (see e.g. equa- 
tion (J3J). Successful prediction of systems with perfect 
ID metal adatom step-edge decoration could facilitate 
the quest for the experimental realization of Luttinger 
liquids in such systems, possibly enabling a convincing 
proof of spin-charge separation with angle resolved pho- 
toemission spectroscopy. 
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